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ABSTRACT 

We present an X-ray analysis of 54 normal elliptical galaxies in the Chandra archive and isolate their 
hot gas component from the contaminating point source emission, allowing us to conduct, for the first 
time, a morphological analysis on the gas alone. A comparison with optical images and photometry 
shows that the hot gas morphology has surprisingly little in common with the shape of the stellar 
distribution. We observe no correlation between optical and X-ray ellipticities in the inner regions 
where stellar mass dominates over dark matter. A shallow correlation would be expected if the gas 
had settled into hydrostatic equilibrium with the gravitational potential. Instead, observed X-ray 
ellipticities exceed optical ellipticities in many cases. We exclude rotation as the dominant factor to 
produce the gas ellipticities. The gas appears disturbed, and hydrostatic equilibrium is the exception 
rather than the rule. Nearly all hydrostatic models can be ruled out at 99% confidence, based of 
their inability to reproduce the optical- X-ray correlation and large X-ray ellipticities. Hydrostatic 
models not excluded are those in which dark matter either dominates over stellar mass inside the 
inner half-light radius or has a prominently cigar-shaped distribution, both of which can be ruled out 
on other grounds. We conclude that, even for rather X-ray faint elliptical galaxies, the gas is at least 
so far out of equilibrium that it does not retain any information about the shape of the potential, and 
that X-ray derived radial mass profiles may be in error by factors of order unity. 



Subject headings: galaxies: cooling flows — galaxies: elliptical and lenticular, cD- 
rays: galaxies — X-rays: ISM 
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1. INTRODUCTION 

The launch of Chandra and XMM-Newton opened 
a new era in our understanding of the hydrodynamic 
histories of galaxies and clusters. Ea rly spectral ev- 
idence from ROSAT and ASC A (e.g. iFabbiano et all 
119941 : iBuote fe Canizaresl Il997t) already suggested soft 
diffuse gas and harder stellar point sources as the two 
dominant components of the X-ray emission of normal 
elliptical galaxies. But Chandra made it possible for the 
first time to spatially resolve a significant fraction of the 
point source component into individual sources, which 
are now believed to c onsist mainly of low - mass X-ray bi- 
naries (LMXBs, e.g. ISarazin et al1l2003l and references 
therein). This has significantly contributed to the un- 
derstanding of the correlation between X-ray and blue 
luminosity. The L-^-Lb diagram shows a steep Lx oc L B 
relation at the gas-dominated group and cluster scale, 
and gets shallower toward low X-ray luminosities for nor- 
mal elliptical galaxies due to th e increasing importanc e 
of the point source component (O'Sulliv an et alj ,2001). 
For these galaxies, the LMXB emission severely contami- 
nates the diffuse hot gas emission and complicates efforts 
to reveal its spatial structure. 

It has long been assumed that the hot interstellar 
medium (ISM) in elliptical galaxies is in hydrostatic equi- 
librium with the un derlying gravitational potential (e.g. 
iForman et al.|[l985l ). The desire to make this assump- 



tion is natural since it then gives us a powerful tool 
to probe the host galaxy's mass distribution. As such, 
radial mass profiles derived from observed X-ray pres- 
sure profiles are among the strongest providers of ev- 
idence for the existence of massive dark matter halos 
surro u nding normal ellipti c al gala x ies fe.g.lForman et alJ 
19851 iKiileen fc Bicknel] Il988t iPaolillo et all T2005 



Humphrey et alj l2006t iFukazawa et al. 2006) and have 
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been found to be consistent wit h those derived from 
gas kinematics in spiral galaxies (jSofue fc Rubinl I2001L 
and refer ences therein) and stellar kinematics in ellipti- 
cals (e.g. IStatler et al.l 11999^ . implying dark halos with 
mass profiles similar to those of isothermal spheres. Un- 
fortunately, stellar kinematics can only probe the dark 
matter content out to a few optical radii. Farther out, 
one has to rely on X -ray mass profiles, gravitational 
lensing (iKeetonl |2001|) or kinemati cal planetary nebula 
(PN) date tiNapolitano et alj|200l) . A recent study by 
iRomanowskv et al.l ( 2003l )~~using PN data advocates a 
rather low dark matter contribution around some ellipti- 
cal g alaxies, at odds wit h standard A-CDM simulations 
(e.g. ISpringel et ai1l2005[ ). This discrepancy may be rec- 
onciled by appeal ing to radial orbits of the halo stars 
(|Dekel et al.ll2005t ). which underlines the difficulty of in- 
terpreting PN kinematics. 

The assumption of hydrostatic equilibrium in elliptical 
galaxies has always been considered well founded and 
understood. That radial mass profiles consistent with 
other methods can be derived from it indicates that it 
is not seriously wrong. If hydrostatic equilibrium holds 
precisely, the isophotes of the hot gas emission should 
very nearly trace the projected potential isophotes, and 
may be used to trace the total mass distribution. Buotc 
and collaborators use this approach to argue for a highly 
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flattened t riaxial dark matter halo in the elliptical galax y 
NGC 720 (|Buote k Canizaresl Il99i IBuote et al.ll200l . 
However, no test has been devised to verify that hy- 
drostatic equilibrium holds precisely enough all the way 
through a galaxy to make this kind of inference valid. 

Indeed, substantial evidence is accumulating to 
the contrary. Combining stellar kinematic informa- 
tion with X-ray o bserva tions of NGC 4472, both 
ICiotti fc Pellegrini! (|2004D and iMathews k Brighentil 
( 2003b!) independently argue for a lack of hy- 
drostatic equilibrium in this particular galaxy. 
IStatler k McNamaral (|2002f ) argue that the extreme 
X-ray flattening of NGC 1700 provides a counterexample 
of an object that cannot be in hydrostatic equilibrium, 
and suggest that it is rotationally flattened. Rotational 
flattening might be expected for a variety of reasons 
(|Mathews k Brighentil l2003al) . A large fraction of the 
hot gas in elliptical galaxies is thought to come from 
stellar mass loss, and should carry the stellar angular 
momentum. A fraction of the hot gas may also be 
acquired externally during mergers, stripped off during 
close enco unters, or fall in from a c ircumgalactic gas 
reservoir (Bri ghenti k Mathews! Il998f ). In each case, 
the gas should contain significant angular momentum. 
In the standard cooling flow scenario the gas should 
slowly flow inward, conserving angular momentum 
and settling into a rot ationall y supp o rted cooling disk 
( Brig henti k Mathews! 119971 Il996t iKlev k Mathews! 



11995ft . However. lHanlan k Bregmanl (|2000f > demon- 
strate, with ROSAT and Einstein data for 6 elliptical 
galaxies, a lack of gas disk signatures. They find 
ellipticities that generally do not exceed values of ~ 0.2, 
whereas disk models predict ellipticities larger than 
-0.5. 

Detailed Chandra and XMM studies of over two dozen 
individual early-type galaxies are now published in the 
literature. Many of these observations reveal systems 
that are morphologically disturbed, with a large va- 
riety of suggested causes. For example, Uones et"al] 
(120021) argue for the presen ce of shocks in NGC 4636; 
iFinoguenov k Jonesl (|2001) find evidence for interac- 
tions with the central radio source in NGC 4374; and 
iMachacek etall (|2004h suggest that NGC 1404 is mov- 
ing relative to the Fornax intracluster gas. Only 
a few obje cts appear round and quie scent, such as 
NGC 4555 (lO'Sullivan k Ponmanl [200l or NGC 6482 
(jKhosroshahi et al.l 120041 )! Nonetheless, the use of 
hydrostatic equilibrium as a fundamental assumption 
about the majority o f early-type galaxies persists (e.g. 
iFukazawa et al.ll2006l ). 

The main obstacle to understanding the true physical 
state of the hot ISM in elliptical galaxies is the separa- 
tion of the gas and unresolved LMXB contributions to 
the diffuse emission, which is particularly important for 
X-ray faint galaxies. So far, two main approaches have 
been employed. (1) Resolved point sources are removed, 
and the residual diffuse image is adaptively smoothed. 
It is then argued or hoped that the contribution of unre- 
solved point sources can be neglected while interpreting 
the diffuse emission a nd that what is sho wn is close to the 
gas mor phology (e.g. IBuote et alJl2002ft . (2) The obser- 
vation is separated into broad radial bins with sufficient 
signal to extract and fit a spectrum with a two compo- 
nent model, thus spectrally separating gas and LMXBs 



in the radial profile (e.g. lHumphrev et al.ll2006f) . Method 
(1) has the obvious disadvantage that one does not know 
what one is looking at. One gains no information on gas 
morphology in X-ray faint systems which are known to be 
LMXB dominated. Method (2) allows some insight into 
the radial distribution and extent of the gas but reveals 
nothing about the true elliptical shape or asymmetric 
features. 

This is the first paper in a series that analyzes data 
on 54 elliptical galaxies in the Chandra public archive. 
We homogeneously reanalyze the observations and intro- 
duce a new technique to isolate the gas from LMXBs . 
An adaptive binning tec hnique (|Diehl k Statlerl 120061 : 
ICappellari k Copinll2003h is then used to reveal the mor- 
phology of the gas alone. We present a gas gallery and a 
quantitative morphological analysis. For the first time, 
we are able to test whether hydrostatic equilibrium gen- 
erally holds in normal ellipticals, by comparing gas and 
stellar isophotal ellipticities in the region where the stel- 
lar mass is expected to dominate the potential. We show 
that the archive sample fails the test. We also exam- 
ine the evidence for rotational support of the gas, and 
find it to be extremely weak. The subsequent papers 
in this series will take up the origin of the observed gas 
morphology (|Diehl k Statlerl l2007al hereafter Paper II) 
and address the importance of ce ntral active galactic nu- 
clei (AGN) in reheating the gas (|Diehl k Statlerl l2007bL 
hereafter Paper III). 

The remainder of this paper is organized as follows. In 
$2]we describe the details of the data reduction pipeline, 
the new isolation technique and the method used to de- 
rive ellipticity profiles. In $3] we present our comparison 
between optical and X-ray properties and address the 
evidence for hydrostatic equilibrium and rotational sup- 
port. We discuss the implications of our findings in [21 
before we summarize in S}5] 

2. DATA 

2.1. Sample Selection and Pipeline Reduction 

We select all galaxies having non-grating ACIS-S ob- 
servations with approved exposure times longer than 
10 ks in the Chandra public archive of cycle 1-4 which 
are classified as E or E/SO in the Lyon-Meud on Extra- 
galactic Database (LEDA: IPaturel et aT1ll997l ). We re- 
move brightest cluster galaxies and objects with purely 
AGN-dominated emission. This sample is almost iden- 
tical to that used in a previous paper (|Diehl k Statlerl 
2005), except for the removal of two galaxies. We ex- 
clude the NGC 4782/NGC 4783 galaxy pair due to its 
ongoing merger and the dwarf elliptical NGC 1705 whose 
luminosity is two orders of magnitude smaller than that 
of the next largest object. Our final sample consists of 
54 early-type galaxies, 34 of which are listed as m embers 
of a g roup in the Lyon Group of Galaxies (LGG. iGarcial 
11993 ) catalog, with 19 identified as the brightest group 
member. Eight are also identified a s members of X-ray- 
brigh t groups in the GEMS survey (jOsmond k Ponmanl 
|2004| ). together with 5 additional galaxies that are not 
listed in the LGG catalog. 

We apply a homogeneous data reduction pipeline to 
all observations starting from their level levent files us- 
ing CIAO version 3.1 with calibration data base 2.28. 
The basic data reduction steps follow the recommenda- 
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tions according to Chandra's ACIS data analysis guide 3 . 
The newest gain file is applied and adjusted for time- 
dependent gain variations to account for the drift of the 
effective detector gains with time caused by changes in 
the charge transfer inefficiency. Observation-specific bad 
columns and pixels are removed and each observation is 
restricted to its good time intervals. We additionally fil- 
ter each light curve by iteratively applying a 2.5cr thresh- 
old to remove background flares. The remaining light 
curve is then clipped at 20% above the average count 
rate, to match the standard of the Markevitch blank sky 
background files 4 . Each light curve is inspected and ver- 
ified to be clean. If the entire observation is affected 
by a very long flare that manifests itself as an under- 
lying "ramp" in the light curve, we flag the object and 
use a local background spectrum in our spectral analysis 
instead of blank sky fields. The background correction 
for the subsequent spatial analysis is not taken from the 
blank sky background files, but rather derived from sur- 
face brightness profile fits (see also W2.2\i which will au- 
tomatically correct for flux offsets from residual flares. 
The remaining event list is filtered to retain standard 
ASCA grades 0, 2, 3, 4, and 6 to optimize the signal- 
to-background ratio. Cosmic ray afterglows, as flagged 
by the CIAO tool acis-detect_afterglow, are removed for 
the purpose of source detection to minimize the number 
of spurious detections 5 . These photons are retained for 
the rest of the analysis. In practice, we find that keeping 
or removing the afterglows makes no difference to the 
results. 

Our analysis is restricted to photon energies between 
0.3-5 keV, maximizing the relative contribution of soft 
hot gas emission, while avoiding the rise of the particle 
background at higher energies. All quoted X-ray lumi- 
nosities are restricted to this band. We split this energy 
range further into a soft (0.3-1.2 keV) and a hard (1.2- 
5.0 keV) band. To create flux-calibrated images that take 
the spatially dependent spectral changes into account, 
we create mono-energetic exposure maps in steps of 7 
in PI (~ 100 eV). The observation is split into 14.6 eV- 
wide (the width of one PI channel) individual images. A 
photon-flux-calibrated "slice" is created by dividing this 
counts image by the energetically closest exposure map. 
The sum of all individual slices represents the final pho- 
ton flux image. These flux-calibrated images allow an 
accurate flux determination even in cases where a spec- 
tral analysis is impossible due to the lack of sufficient 
signal. 

2.2. Isolating the X-ray Gas Emission 

To isolate the hot gas emissio n we f ollow the proce- 
dure outlined bv lDiehl & Statleri (|2005l ). A summary is 
repeated here with some additional details, which are 
essential for the morphological analysis in this paper. 

We use the CIAO tool wavdetect to identify point 
sources, remove regions enclosing 95% of the source flux 
and refill the holes with simulated Poisson counts to ob- 
tain an image of diffuse emission. We determine uniform 
background values for the soft and hard bands by ex- 
tracting radial surface brightness profiles for the diffuse 

http: / / cxc.harvard.edu / ciao / guides / acis_data.html 



emission from the calibrated photon-flux images and fit- 
ting them with /?, double-/? and Sersic models plus an 
additive constant background. The sole purpose of these 
fits is to determine the background value; all other fitting 
parameters are discarded for the subsequent analysis. As 
the models generally yield slightly different background 
values, we add half of this difference to the statistical 
background error. If there is insufficient signal to pro- 
duce an accurate fit, we compute the average surface 
brightness level outside a 2.5 arcmin radius and use it 
as the background value. We then subtract the spatially 
uniform backgrounds from the photon-flux calibrated im- 
ages to obtain the background corrected images of diffuse 
emission for the soft and hard bands. Background er- 
rors are propagated into the final gas variance image and 
taken into account in the subsequent analysis. 

The hot gas emission is still contaminated by unre- 
solved point sources and the incompletely removed PSF 
wings of resolved point sources. These contaminants 
must be modeled and removed to isolate the hot gas emis- 
sion alone. To do this, we split the background corrected, 
photon-flux calibrated images into the soft band S and 
hard band H. Gas and point source components con- 
tribute at different levels to each band, determined by 
their respective softness ratios 7 and S. By expressing 
both bands as linear combinations of the hot gas com- 
ponent G and the unresolved point source component P, 
we can solve this system of equations to isolate the gas 
emission itself: 



thus, G 



S = jP + SG, 
ff=(l- 7 )P 

1-7 



(5 — 7 



S 



1-5)G 

7 



1 



7 



H 



(1) 
(2) 

(3) 



This decomposition depends on the softness ratios 7 and 
5. To determine the softness of the unresolved point 
sources, we first analyze the hardness ratios of resolved 
point sources. We find no evidence for any spatial depen- 
dence or luminosity dependence of their spec tral proper- 
ties. This is in agreement with studies by llrwin et alJ 
( 2003]), which suggest a universal nature of LMXBs, and 
allows us to use the known spectral properties of resolved 
LMXBs as a template for their unresolved counterparts. 
For each galaxy, we select resolved point sources between 
5" and 5 optical effective radii from the center to avoid 
the influence of the central AGN, and to minimize con- 
tributions from serendipitous background sources. We 
exclude high luminosity sources (> 200 counts) to en- 
sure that the spectral fits are driven by low-luminosity 
sources. If more than 10 sources fulfill these selection 
criteria, we fit an absorbed power law model to the com- 
bined LMXB spectrum. We find that a narrow range 
of photon indices (1.1 to 1.8) is adequate for the point 
source component in all of our galaxies) . To get an equiv- 
alent point source spectral model for galaxies without 
sufficient resolved point sources, we simultaneously fit a 
power-law to all low-luminosity (Lx < 5 x 10 37 erg s~ x ) 
LMXBs available in our complete galaxy sample. This 
combined fit yields a photon index of 1.603, which we 
adopt as the representative spectral model for source- 
poor galaxies. By integrating the adopted unabsorbed 



4 http://cxc.harvard.edu/cal/ Acis/CaLprods/bkgrnd/acisbg/COOKB?R9I^ral model over the soft and hard band, we derive 

5 http://cxc.harvard.edu/ciao/threads/acisdetectafterglow/ the point Source softness ratio 7. 
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To find the softness ratio of the hot gas component, S, 
we extract a spectrum from the diffuse emission within 
only the inner 3 optical effective radii of each galaxy, to 
avoid contamination by intragroup or intracluster gas. 
We use an absorbed single temperature APEC thermal 
plasma model for the hot gas and add the adopted power- 
law model to represent unresolved point sources. We fix 
the redshift parameter to the LEDA value, while tem- 
perature, metallicity and normalizations are generally 
allowed to vary freely. For low signal-to-noise spectra, 
we fix the metallicity to the solar value. All described 
spectral fits are computed with the CIAO tool Sherpa 
and corrected for Galactic absorption, with the hydrogen 
column density fixed at the Galactic value for the line of 
sight, as determined by the CIAO tool Colden 6 . This 
decomposition technique removes both the unresolved 
point sources and the PSF wings of resolved sources, due 
to their universal spectral nature. 

To reliably interpret the sparse gas images and to re- 
veal any spatial features, we bin the gas images with an 
adapti ve binning method using weighted Voronoi tessela - 
tions (iDiehl fe Staffer! 12006: Ca ppellari fe Cot)inl[2003h 7 
to achieve an approximately constant signal-to-noise ra- 
tio per bin of 4. In background dominated, very-low 
signal-to-noise regions, the required bin size can get very 
large and the bins occasionally "eat" into the inner emis- 
sion to accrete sufficient signal. To avoid this effect, we 
restrict the maximum bin size for these objects, resulting 
in a drop of signal-to-noise per bin in the outer regions. 

A few things should be kept in mind regarding our 
new isolation technique for the hot gas. First, we as- 
sume isothermal gas throughout the galaxy, resulting in 
a spatially constant 5 parameter. We test the validity 
of this assumption by creating two-dimensional temper- 
ature maps for systems with the highest signal-to-noise 
data and correct the gas image for a spatially dependent 
S value. Although most galaxies exhibit temperature gra- 
dients, the corrections to the gas surface brightness are 
generally < 10% and do not affect the overall gas mor- 
phology. For those objects that do have sufficient sig- 
nal to produce a spatially dependent temperature map, 
the effects on ellipticity and asymmetry are negligible. 
Spatial temperature gradients and inhomogeneities will 
be examined in more detail in Paper III. However, one 
should be aware that strong localized temperature dif- 
ferences can potentially result in an over-subtraction of 
hotter features and under-subtraction of colder features. 

Second, the region influenced by the central PSF of an 
AGN may be subtracted incorrectly, as the AGN softness 
ratio is most likely different from that of the unresolved 
point source component. Thus, we exclude the central 
regions from our analysis, and one should refrain from 
interpreting the gas maps at the very center. A detailed 
analysis of the central AGN and its effects on the overall 
morphology will be presented separately in Paper II. One 
should also keep in mind that possible direct X-ray jet 
signatures could be masked out in our analysis, as jet 
knots would most likely be identified as point sources 
and thus removed. 8 

s http:/ /cxc. harvard.edu/toolkit/colden.jsp 

7 http://www.phy.ohiou.edu/~diehl/WVT 

8 NGC 315 is the only object in our sample with a conspicu- 
ous jet. We have verified by eye that the jet is removed by our 



A drawback of our approach is the inevitable loss of 
signal. Since we are subtracting a scaled version of the 
hard band from the soft band to isolate the gas emis- 
sion, we effectively subtract a part of the gas flux. Even 
though our technique corrects for this missing flux, one 
still loses a certain amount of signal, reducing the spatial 
resolution in the adaptively binned gas image. The soft- 
ness of the gas comes close to that of the unresolved point 
sources for temperatures exceeding ~ 1.5 keV. From that 
point on, both components become increasingly indistin- 
guishable and the signal-to-noise ratio of the gas image 
can be very low, even if there is substantial gas present. 
Fortunately, few objects in our sample have temperatures 
exceeding 1 keV. 

The definite advantage of our algorithm is that we can 
properly isolate the hot gas emission, without assump- 
tions about the spatial distribution of unresolved point 
sources, or their total flux. Even the effect of a spatially 
variable point source detection limit (due to variations 
in background levels and the changes in size and shape 
of the point spread function) is compensated for by our 
algorithm. 

To give the reader a better understanding of the rela- 
tive contributions of the hot gas and the resolved and 
unresolved point source components, we list the rela- 
tive fractions in Table Q] for a radial range between 2.5" 
and 3 J-band effective radii This shows that for a 

large subset of galaxies, a significant percentage of point 
sources is still unresolved and contaminates the diffuse 
emission, disguising the true gas morphology. In some 
cases (e.g. NGC 3115) the diffuse emission is even con- 
sistent with being almost entirely due to unresolved point 
sources. 

2.3. X-ray Gas Luminosity 

To calculate the total X-ray gas luminosity, we bin the 
calibrated gas images into adaptively sized circular an- 
nuli and p roduce radial surface b rightness profiles, as de- 
scribed in IDiehl fe Statlen (|2005l ). X-ray surface bright- 
ness profiles of galaxy clusters and groups are generally 
well described by (3 or double-/3 profiles with (3 values be- 
tween ~ 0.5 — 1.0. Our (3 model fits indicate that the gas 
profiles of normal ellipticals are generally shallower. For 
28 out of 45 galaxies with sufficient signal to constrain 
the fit, the best fit (3 values are < 0.5, consistent with 
the known relation betw een gas tempera ture and (3 for 
groups and clusters (e.g. IVoit et~ai1l2002D . 

However, a (3 value below 0.5 yields infinite total flux 
when extrapolated to large radii. As this renders the 
(3 models unusable for determining total luminosities, 
we adopt the well-known Sersic models instead. The 
Sersic model has the advantage of yielding finite fluxes 
and produces equivalently good fits. The values for our 
best Sersic model fits are published in an earlier paper 
(jDiehl fc Statlerll2005h . We derive X -ray gas luminosi- 
ties by summing the calibrated gas images over the field 
of view, and use the model fits to correct for missing 
flux outside the field of view. The luminosities listed in 
Table[T]are also corrected for absorption effects by multi- 
plying a correction factor, derived from the best spectral 
fit, which we integrate over energy with and without ab- 
sorption by the Galactic neutral hydrogen. The error 

procedure. 
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TABLE 1 
Chandra X-ray properties 



N £11116 


Obsld 


r a 


fr b 

J Gas 


f b 

J unrcs 


/resol 


b 








PA Y d 


IC1262 


2018 


31 


58 ± 15% 


41 ± 15% 


i ± 


1% 


2.0 ± 1.7 x 


10 43 


0.43 ± 0.09 


69 ± 11 


IC1459 


2196 


52 


58 ± 10% 


12 ± 10% 


31 ± 


1% 


4.3 ±3.2 x 


1Q 40 


0.16 ±0.08 


32 ± 20 


IC4296 


3394 


24 


74 ± 10% 


7± 9% 


18 ± 


1% 


1.1 ± 0.4 x 


10 41 


0.28 ±0.06 


55 ± 14 


NGC0193 


4053 


29 


85 ± 13% 


10 ± 9% 


5± 


1% 


2.5 ±0.8 x 


10 41 






NGC0315 


4156 


52 


58 ± 10% 


32 ± 10% 


10 ± 


1% 


9.4 ±3.4 x 


10 40 


0.13 ±0.16 


14 ±32 


NGC0383 


2147 


43 


49 ± 10% 


35 ± 10% 


16 ± 


1% 


< 7.5 x 


10 41 


0.23 ±0.09 


114 ±25 


NGC0404 


870 


24 


94 ± 43% 


± 42% 


6± 


2% 


< 2.1 x 


10 38 






NGC0507 


317 


26 


69 ± 13% 


30 ± 13% 


1 ± 


1% 


> 5.7 x 


10 42 


0.13 ±0.17 


91 ± 33 


NGC0533 


2880 


35 


89 ± 19% 


9 ± 18% 


2± 


1% 


9.6 ±3.5 x 


10 41 


0.37 ± 0.05 


26 ±3 


NGC0720 


492 


34 


74 ± 10% 


15 ± 9% 


12 ± 


1% 


9.3 ±2.7 x 


10 40 


0.09 ± 0.07 


75 ± 28 


NGC0741 


2223 


29 


58 ± 11% 


21 ± 11% 


22 ± 


1% 


3.2 ± 1.3 x 


10 41 


0.24 ± 0.05 


162 ± 13 


NGC0821 


4006 


13 


17 ± 23% 


5 ± 87% 


78 ± 54% 


< 3.3 x 


10 40 






NGC1132 


801 


12 


83 ± 16% 


16 ± 16% 


1 ± 


1% 


> 9.1 x 


10 42 


0.08 ±0.24 


40 ±46 


NGC1265 


3237 


82 


15 ± 8% 


38 ± 8% 


47 ± 


5% 


< 1.1 X 


10 42 






NGC1316 


2022 


26 


78 ± 13% 


5 ± 13% 


18 ± 


1% 


5.7 ± 2.1 x 


10 40 


0.33 ± 0.06 


14 ± 13 


NGC1399 


319 


24 


85 ± 23% 


7 ± 23% 


8± 


1% 


> 7.9 x 


10 41 


0.34 ± 0.04 


179 ± 7 


NGC1404 


2942 


29 


94 ± 16% 


1 ± 16% 


5± 


1% 


1.7 ±0.4 x 


10 41 


0.06 ±0.07 


63 ± 26 


NGC1407 


791 


43 


70 ± 10% 


18 ± 10% 


12 ± 


1% 


1.0 ±0.3 x 


10 41 


0.13 ±0.06 


165 ± 16 


NGC1549 


2077 


20 


52 ± 10% 


18 ± 9% 


31 ± 


2% 


> 2.0 x 


10 40 






NGC1553 


783 


19 


61 ± 13% 


24 ± 13% 


15 ± 


1% 


2.8 ± 2.6 x 


10 40 


0.57 ±0.14 


90 ± 15 


NGC1600 


4283/4371 


49 


79 ± 13% 


17 ± 12% 


4± 


1% 


> 1.2 x 


10 42 


0.24 ± 0.09 


178 ± 17 


NGC1700 


2069 


39 


85 ± 14% 


7 ± 14% 


9± 


1% 


> 3.2 x 


10 41 


0.20 ±0.06 


89 ± 12 


NGC2434 


2923 


25 


69 ± 9% 


9± 6% 


21 ± 


1% 


2.6 ±2.0 x 


10 40 






NGC2865 


2020 


25 


34 ± 21% 


60 ± 19% 


6± 


2% 


< 9.9 x 


10 40 






NGC3115 


2040 


14 


2± 9% 


35 ± 8% 


63 ± 


5% 


< 8.7 x 


10 39 






NGC3377 


2934 


39 


15 ± 22% 


± 15% 


85 ± 10% 


< 6.1 x 


10 39 






NGC3379 


1587 


30 


14 ± 4% 


11 ± 4% 


75 ± 


2% 


< 6.3 x 


10 39 






NGC3585 


2078 


35 


32 ± 7% 


19 ± 5% 


48 ± 


3% 


> 4.2 x 


10 39 






NGC3923 


1563 


16 


74 ± 9% 


11 ± 8% 


15 ± 


1% 


4.3 ± 1.3 x 


1Q 40 


0.15 ±0.10 


44 ±25 


NGC4125 


2071 


63 


80 ± 13% 


1 ± 13% 


19 ± 


1% 


7.2 ± 2.7 x 


1Q 40 


0.24 ±0.11 


71 ± 10 


NGC4261 


834 


31 


63 ± 11% 


18 ± 10% 


19 ± 


1% 


4.8 ± 1.1 x 


1Q 40 


0.08 ±0.08 


108 ± 35 


NGC4365 


2015 


40 


44 ± 8% 


14 ± 8% 


41 ± 


1% 


> 3.8 x 


10 40 






NGC4374 


803 


28 


86 ± 18% 


5 ± 18% 


9± 


1% 


5.9 ± 1.3 x 


1Q 40 


0.38 ±0.06 


116 ± 7 


NGC4406 


318 


14 


89 ± 12% 


10 ± 12% 


2± 


1% 


> 1.0 x 


10 42 


0.17 ±0.10 


113 ± 24 


NGC4472 


321 


35 


83 ± 16% 


12 ± 16% 


5± 


1% 


> 8.5 x 


10 41 


0.12 ±0.04 


136 ± 13 


NGC4494 


2079 


19 


0± 2% 


43 ± 10% 


57 ± 


5% 


< 2.1 x 


10 40 






NGC4526 


3925 


38 


55 ± 12% 


8 ± 10% 


37 ± 


2% 


8.8 ± 7.5 x 


10 39 


0.07 ±0.15 


1 ±43 


NGC4552 


2072 


53 


80 ± 10% 


4 ± 10% 


16 ± 


1% 


2.1 ± 1.2 x 


1Q 40 


0.21 ±0.03 


7±4 


NGC4555 


2884 


27 


73 ± 20% 


19 ± 18% 


9± 


2% 


> 2.3 x 


10 41 






NGC4564 


4008 


17 


72 ± 23% 


17 ± 15% 


11 ± 


2% 


> 2.0 x 


10 39 






NGC4621 


2068 


24 


42 ± 11% 


7± 9% 


50 ± 


4% 


1.1 ±0.9 x 


1Q 40 






NGC4636 


323 


46 


95 ± 17% 


3 ± 17% 


3± 


1% 


2.7 ±2.0 x 


10 41 


0.35 ±0.04 


9±3 


NGC4649 


785 


31 


80 ± 11% 


11 ± 11% 


9± 


1% 


1.3 ±0.3 x 


10 41 


0.08 ±0.03 


95 ± 27 


NGC4697 


784/4727/4728 


112 


57 ± 6% 


6± 5% 


38 ± 


1% 


> 3.5 x 


10 40 


0.20 ±0.16 


80 ±23 


NGC5018 


2070 


25 


61 ± 11% 


21 ± 9% 


18 ± 


2% 


< 1.9 x 


10 41 






NGC5044 


798 


20 


100 ± 18% 


± 18% 


0± 


1% 


2.6 ±0.8 x 


10 42 


0.41 ±0.08 


27 ± 10 


NGC5102 


2949 


34 


50 ± 22% 


4 ± 13% 


46 ± 


7% 


< 1.6 x 


10 39 






NGC5171 


3216 


34 


83 ± 25% 


13 ± 20% 


4± 


1% 


> 2.7 x 


10 42 






NGC5532 


3968 


48 


28 ± 6% 


24 ± 5% 


48 ± 


2% 


< 8.7 x 


10 41 






NGC5845 


4009 


30 


57 ± 30% 


41 ± 19% 


2± 


3% 


< 5.2 x 


10 40 






NGC5846 


788 


24 


95 ± 16% 


4 ± 16% 


2 ± 


1% 


3.9 ±0.9 x 


10 41 


0.22 ±0.05 


175 ± 8 


NGC6482 


3218 


19 


90 ± 13% 


9 ± 13% 


1 ± 


1% 


1.7 ± 1.3 x 


10 42 


0.18 ±0.09 


36 ± 23 


NGC7052 


2931 


9 


90 ± 14% 


6 ± 10% 


4± 


1% 


> 1.1 X 


10 41 


0.26 ±0.12 


54 ±23 


NGC7618 


802 


11 


85 ± 11% 


14 ±11% 


1 ± 


1% 


2.3 ±0.9 x 


10 42 


0.38 ±0.14 


99 ±16 



a Effective Chandra exposure time in ks after removal of flares. b Percentages of observed photon fluxes for the 
hot gas, unresolved and resolved point sources, integrated between 2.5" and 3Rj . c Total X-ray gas luminosity in 
ergss -1 for the 0.3 — 5keV band. d X-ray gas ellipticity and position angle, evaluated between 0.6 and 0.9Rj. 
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bars for luminosities can get quite large for objects with 
very wide gas emission, where a significant fraction of 
flux is derived from the uncertain extrapolation to large 
radii. The quoted uncertainties in Table [1] also include 
systematic errors associated with uncertainties in the 5 
and 7 parameters from section 12.21 as well as the un- 
certainty in our adop ted distances (for more details, see 
iDiehl fe StatTeHl2005l) . In cases where the uncertainties 
are larger than the actual values due to uncertainties in 
the extrapolation or where we have insufficient signal for 
a surface brightness model fit, we determine the lumi- 
nosity by summing the flux in the field of view. If the 
summed flux represents a detection at a > 3tr level, we 
report the 3<r lower bound as a low limit in Table [TJ If 
this is not the case, we report the 3<r upper bound as an 
upper limit. 

2.4. X-ray Ellipticity and Position Angle Profiles 

The extremely sparse nature of X-ray data generally 
prohibits the use of isophote fitting techniques commonly 
employed with optical data. Moreover, because our gas- 
only images are obtained from scaled differences of the 
hard and soft bands, they contain many individual pix- 
els with large negative counts, which can render standard 
algorithms numerically unstable. To avoid this problem, 
we fit isophotes to the adaptively binned gas images, 
adapting techniques used in TV-body simulations. 

We populate each bin randomly with a number of 
pseudo counts (particles) such that the expected >/N 
Poisson fluctuations match the signal-to-noise ratio in 
the bin. The pseudo count fluxes (masses) are chosen 
to give the correct flux in each bin. Isophotes are then 
fitted to the pseudo count distribution using an iterative 
algorithm that diagonalizes the second-moment tensor in 
a thin elliptical ring and manipulates the axis ratio until 
the eigenvalues ma tch those for a constant de nsity ring 
of the same shape (jStatler fc McNamarall2002l ). Twenty 
random realizations are run, and the distribution of com- 
plex ellipticities is used to find the mean isophotal el- 
lipticity, major axis position angle, and errors at each 
radius. 

Tests show that this technique is able to robustly re- 
cover simulated isophotal profiles. An example is shown 
in figure [TJ Here we have simulated soft- and hard-band 
images from separate gas and point source contributions. 
Both gas and point sources follow different radial sur- 
face brightness models and have very different isopho- 
tal shapes. The point sources are distributed as a de- 
Vaucouleurs profile with an effective radius of 50 pixels 
with a constant ellipticity of 0.1 and major-axis position 
angle of —10°. We draw the point source luminosities 
from a typical power-law luminosity function with a log- 
arithmic slope of —2.0. We then compute the Chandra 
point spread function with the CIAO tool "mkpsf" at 
the source location and finally simulate Poisson counts. 
The gas on the other hand follows a (3 profile with a core 
radius of 20 pixels and (3 = 2/3. The more complicated 
gas ellipticity and position angle profiles are shown by 
the smooth curves in the middle and bottom panels of 
Figure [TJ We assume typical softness ratios of 7 = 0.5 
and 5 = 0.9 and choose the luminosities such that 60% 
of the counts in the simulated full-band image are from 
gas. This percentage is representative of the lowest end 
of the range for the objects in Table I for which isophote 




Semimajor axis (pixels) 

Fig. 1. — Tests of isophote fitting for simulated data with 40, 000 
(red), 10, 000 (green), and 2, 500 (blue) counts from gas. Gas emis- 
sion is separated from point source emission by identifying and re- 
moving resolved sources and differencing soft- and hard-band im- 
ages to remove unresolved sources (§ 2.2) and PSF wings. Isophotes 
are fitted to the adaptively binned gas image as described in § 2.4. 
Top panel shows gas surface brightness (in counts/pixel). Mid- 
dle and bottom panels show isophotal ellipticity and major axis 
poisition angle. Smooth curves indicate the input profiles. 

fitting is possible. A flat background of 0.01 counts/pixel 
is added to the soft and hard images before simulating 
discrete counts. We then go through our entire analy- 
sis pipeline: detecting point sources, removing resolved 
point sources, differencing the images to extract the gas 
map, adaptively bin, and fit isophotes as described above. 
Figure [TJ shows the results for 40, 000, 10, 000, and 2, 500 
gas counts for a comparable radial range as is used for the 
real data. The isophotal profiles are recovered to within 
the statistical errors. Clearly, at least a few thousand 
counts in gas are necessary to recover accurate elliptic- 
ities. Only in cases where the entire central region is 
identified as consisting of several point sources and is 
completely blocked out, we note a systematic deviation 
from the input profiles. In our data this region is found 
to be within the central 10 pixels. Thus, we do not inter- 
pret the central 5" (10 pixels) in our isophotal profiles. 
Furthermore, accurate isophotes can be fit only where 
the gas surface brightness is > 2 times the background. 
As a result, the isophotes are completely insensitive to 
errors in the background level. 

2.5. Optical Data 

Optical data for the sample galaxies are given in Table 
[2j We adopt the effective J-band ellipticity ej, position 
angle PAj and half-light radius Rj (Table ^ , as well as 
the absolute K magnitude from t he 2MASS extended 
source catalog ([Jarrett et al.l l2000h . These ellipticities 
and position angles are extracted at a constant mag- 
nitude level, approximately 3cr above background noise 
level, generally translating to an extraction radius be- 
tween 20-40 arcsec. 

In addition, we extract i?-band images from the sec- 
ond epoch of the Digitized Sky Survey (DSS-2R) and 
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make use of publicly available optical surface photom- 
etry from the literature. In particular, we adopt B 
(20 galaxies), V (20) and I (19) photometry from 
iGoudfrooii etail (fl99l 9 V (21), R (25), and / (23) 

U (13), B (15), and R (14) 
and F 814VF (9) photometry 



Bender et all (|198. 
Peletier ct al. ( 199 



Falcon-Barroso et al. (2006). Refer to Tabled for 



from 
from 
from 

details. We use the optical photometry to extract ef- 
fective optical position angles PA opt and ellipticities e op t 
between 0.6-0.9 Rj. 

We also query the Hyperleda (jPrugniel et alj 119981 ) 
data base to derive average rotational velocities between 
0.6-0.9 R j. In cases where Hyperleda does not provide 
the data in electronic form, we read off approximate val- 
ues from published kinematic profile plots. The list of 
references that are used to derive the adopted rotational 
velocity is given in Table [2] This rotational velocity is 
closely related to the maximal rotational velocity but is 
generally a better defined quantity, as it is independent of 
the observational cutoff of the available kinematic data. 

3. RESULTS 
3.1. X-ray vs. Optical Luminosity 

For the first time, we have spatially isolated the hot gas 
emission in normal elliptical galaxies to allow an unbi- 
ased morphological analysis of a large sample of objects. 
Figure [5] shows the X-ray gas gallery, along with the op- 
tical DSS-2 i?-band images for each object. We rescale 
all images to show the same physical size of 50 kpc, ex- 
cept for the two closest galaxies, where we indicate the 
size by a separate scale bar. We have ordered Figure [5] 
with decreasing absolute K luminosity, a good indicator 
of the total stellar mass of the host galaxies. This is an 
intuitive way to visualize t he large scatter in the well- 
known Lx oc L| relation (O'Sulliv an et al.l I2001L and 
references therein), one of the elementary empirical cor- 
relations connecting X-ray and optical parameters. Its 
scatter shows up as a large spread in X-ray properties 
for galaxies that have similar stellar contents, i.e. that 
are close to each other in the figure. Galaxies with simi- 
lar optical properties can have gas luminosities differing 
by orders of magnitude (e.g. NGC 5044 and NGC 1549, 
NGC 4555 and NGC 0533). At the same time, galaxies 
with comparable gas contents can have extremely differ- 
e nt optical appearances ( e.g. NGC 1404 and NGC 4649). 

lO'Sullivan et al.l (2001) also observe a flattening of the 
Lx~Lb relation toward lower blue luminosities. They at- 
tribute this trend to the increasing relative importance 
of point source emission, which is unresolved in their 
ROSAT and Einstein data. Our Lx,Gas~^B relation, 
with the stellar source contribution removed (Figure |4|), 
is consistent with a si ngle slope, and agrees well with 
lO'Sullivan et all (|2001h 's best fit (dashed line). Nonethe- 
less, we still see a large scatter in this relation, cover- 
ing almost two orders of magnitude in X-ray luminos- 
ity foi_a_gjven_jbsolute_^B magnitude. A recent study 
bv lEUis &: O' Sullivan! (|2006f ) shows a similar relation be- 
tween X-ray and i-T-band luminosities, with an almost 
identical scatter, suggesting a large scatter in X-ray lu- 
minosities as the dominant cause. 



3.2. X-ray vs. Optical Morphology 

Not only do the luminosities of the hot gas and stellar 
components differ vastly, there are also significant differ- 
ences between the X-ray and optical morphology. Figure 
[5] shows that it is impossible to predict the gas morphol- 
ogy from looking at the stellar distribution, or vice versa. 
There are some optically flat galaxies with round X-ray 
isophotes (e.g. NGC 0720), and others with very flat X- 
ray emission (e.g. NGC 1700). Optically round galaxies 
show a similar range in X-ray ellipticities, ranging from 
round (e.g. NGC 1404) to flat (e.g. NGC 0533). 

Many galaxies, in fact, appear asymmetrically dis- 
turbed. We discuss this asymmetry and its origin in 
Paper II. For now, it is important to first quantify 
the overall shape of the gas emission. If the gas is 
in true hydrostatic equilibrium, the hot gas isophotes 

should v ery nearly trace t he pr ojected equipoten- 

tials (e.g. iBinnev fc Tremainelll987t iBuote fc Canizaresl 
|1994|) . Thus, we construct X-ray gas ellipticity and posi- 
tion angle profiles, shown in Figure [5] (solid circles), with 
published optical profiles overlaid. The optical profiles 
are generally "well-behaved" and reveal only modest ra- 
dial trends in their ellipticities and position angles. Op- 
tical isophotal twists, if present at all, rarely exceed a 
few degrees. 

The X-ray profiles tell a completely different story. Of- 
ten, the gas ellipticities change rapidly as a function of 
radius (e.g. NGC 1399, NGC 4374). These features 
are often accompanied by sudden changes in position 
angle (e.g. NGC 4636, NGC 5044). Isophotal twists 
in the gas emission are common, and only very few ob- 
jects show profiles that are consistent with a single ma- 
jor axis orientation (e.g. NGC 3923, NGC 1700). In 
many cases isophotal twists or sudden changes in ellip- 
ticity also coincide with asymmetries in the gas images 
(e.g. NGC 5044). 

It is generally agreed that the potentials of early- 
type galaxies are stellar-mass d ominated inside 1 or 
2 op tical effective radii (e.g . lLoewenstein fe White! 
19991 iMamon fe Lokasl 120051 : iLintott et al.l 120061 : 



Humphrey fc Buotd 120061 ) . To test whether the X-ray 



9 PA profiles from Goudfrooii ct al. (1994) have been corrected 
for mirror-image flips by matching the profiles to other published 
photometry and DSS-2R images. 



gas flattening is consistent with hydrostatic equilibrium, 
we extract mean ellipticities for the hot gas emission and 
the starlight between 0.6 and 0.9 Figure [3] shows 
expanded views of the central regions of the galaxies 
for which gas ellipticities are measurable in this region, 
compared with 2MASS images. The radial range is 
chosen to maximize the available number of objects with 
valid optical and X-ray ellipticity profiles in the same 
interval, resulting in a subset of 24 galaxies. This subset 
spans more than an order of magnitude in optical lumi- 
nosity and nearly two orders of magnitude in X-ray gas 
luminosity. Optical ellipticities are representative of the 
general population, and X-ray gas fractions range from 
57% to 100%. Figure [6] shows that there is absolutely 
no correlation between optical and gas ellipticities. We 
have verified that this result holds for other choices of 
extraction annuli ranging between 0.5 and 1.2 and 
for different choices of optical bandpass. The large scat- 
ter is not caused by systems with obviously disturbed 
X-ray morphology. Of the 8 most highly disturbed 
systems (quantified by the asymmetry index defined 
in Paper II), 6 do not have published optical surface 
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TABLE 2 
Optical properties 



2MASS b LEDA° Hypcrlcda d Photometry" 



Name 


D 






M K 








<=J 


PAj 


Lb 






v Iot Ref. 




Copt 




PA opt 






Ref. 


IC1262 


143.8 ± 


21.6 


-25 


.43 ± 


.33 


14. 


4 


0.34 


85 


5.4 ± 1.7 


X 


10 10 




















IC1459 


29.2 ± 


3.8 


-25 


.53 ± 


.28 


29. 


.1 


0.19 


42 


4.7 ± 2.3 


X 


10 10 






0. 


27 ± 


01 


36.6 ± 


0. 


7 


T 


IC4296 


51.6 ± 


7.7 


— 26. 


.06 ± 


.33 


25. 


5 


0.06 


70 


1.3 ± 0.4 


x 


10 11 






0. 


13 ± 0. 


01 


65.9 ± 


1. 


{) 


u 


NGC0193 


60.8 ± 


9.1 


-24 


.71 ± 


.33 


14. 


5 


0.20 


55 


3.0 ± 1.0 


x 


10 10 




















NGC0315 


71.9 ± 


10.8 


-26 


.33 ± 


.33 


22. 





0.22 


50 


1.4 ± 0.4 


x 


10 11 


Ill 


,\ 


0. 


27 ± 


01 


42.5 ± 


0. 


5 


V 


NGC0383 


73.2 ± 


11.0 


-25 


.84 ± 


.33 


17. 


.8 


0.16 


15 


1.0 ± 0.5 


x 


10 11 


75 


A 
















NGC0404 


3.3 ± 


0.2 
















7.0 ± 1.0 


x 


10 8 




















NGC0507 


71.7 ± 


10.8 


-25 


.98 ± 


.33 


20. 


.1 


0.12 


75 


1.1 ± 0.5 


X 


10 11 




















NGC0533 


77.6 ± 


11.6 


-26 


.01 ± 


.33 


25. 


2 


0.20 


00 


1.0 ± 0.5 


X 


10 11 


20 


B 
















NGC0720 


27.7 ± 


2.2 


-24 


.94 ± 


.17 


27. 


4 


0.41 


145 


3.6 ± 1.4 


X 


10 10 






0. 


41 ± 


01 


140.7 ± 


0. 


1 


I' 


NGC0741 


79.0 ± 


11.8 


-26 


.19 ± 


.33 


25. 





0.38 


95 


1.6 ± 0.5 


X 


10 11 


13.3 


Q 


0. 


18 ± 0. 


01 


87.8 ± 


1. 


2 


V 


NGC0821 


24.1 ± 


1.9 


-24 


.01 ± 


.17 


23. 





0.24 


25 


2.5 ± 0.9 


X 


10 10 


09 


rj 


0. 


38 ± 0. 


01 


32.3 ± 


0. 


fi 


w 


NGC1132 


98.2 ± 


14.7 


-25 


.70 ± 


.33 


19. 


8 


0.30 


140 


8.4 ± 4.2 


X 


10 10 




















NGC1265 


109.5 ± 


16.4 
















2.1 ± 0.9 


X 


10 11 




















NGC1316 


21.5 ± 


1.7 


-26 


.07 ± 


.17 


49. 


8 


0.30 


52 


9.8 ± 9.2 


X 


10 10 


134 


















NGC1399 


20.0 ± 


1.5 


-25 


.19 ± 


.16 


36. 





0.04 


150 


4.9 ± 1.6 


X 


10 10 


31 


Ji 1 


0. 


10 ± 0. 


01 


107.4 ± 


1. 





I' 


NGC1404 


21.0 ± 


1.8 


-24 


.79 ± 


.19 


19. 


3 


0.12 


163 


3.2 ± 0.7 


X 


10 10 


89 


p 


0. 


14 ± 


.01 


160.1 ± 


1. 


2 


T 


NGC1407 


28.8 ± 


3.5 


-25 


.60 ± 


.26 


36. 


4 


0.07 


20 


7.6 ± 3.6 


X 


10 10 


< 20 


Q 


0. 


04 ± 


01 


51.8 ± 


0. 


5 


I' 


NGC1549 


19.7 ± 


1.6 


-24 


.69 ± 


.18 


29. 





0.10 


143 


3.4 ± 0.7 


X 


10 10 


GO 


Q 


0. 


13 ± 


.01 


119.2 ± 


2. 


5 


I' 


NGC1553 


18.5 ± 


1.5 


-25 


.06 ± 


.17 


33. 





0.34 


158 


4.4 ±1.1 


X 


10 10 


188 


JJ 
















NGC1600 


66.0 ± 


9.0 


-26 


.06 ± 


.33 


24. 


8 


0.28 


5 


1.3 ± 0.5 


X 


10 11 


< 10 


rj 


0. 


34 ± 


.01 


7.4 ± 


0. 


6 


V 


NGC1700 


54.4 ± 


8.2 


-25 


.59 ± 


.33 


15. 





0.30 


90 


8.2 ± 3.2 


X 


10 10 


03 


fj 


0. 


26 ± 0. 


01 


88.1 ± 


0. 


5 


I' 


NGC2434 


21.6 ± 


2.9 


-23 


.78 ± 


.29 


19. 


3 


0.08 


145 


2.2 ± 0.6 


X 


10 10 


10 


I 
















NGC2865 


37.8 ± 


3.5 


-24 


.43 ± 


.20 


14. 


8 


0.22 


155 


3.4 ± 1.0 


X 


10 10 


85 


,J 
















NGC3115 


9.7 ± 


0.4 


-24 


.05 ± 


.09 


36. 


4 


0.57 


42 


1.6 ± 0.9 


x 


10 10 


254 


















NGC3377 


11.2 ± 


0.5 


-22 


.81 ± 


.09 


27. 


7 


0.40 


42 


8.0 ± 1.5 


x 


10 9 


80 




0. 


50 ± 0. 


01 


41.0 ± 


0. 


2 


w 


NGC3379 


10.6 ± 


0.5 


-23 


.85 ± 


.11 


29. 





0.09 


70 


1.6 ± 0.3 


x 


10 10 


SO 




0. 


11 ± 


01 


70.6 ± 


1. 


2 


w 


NGC3585 


20.0 ± 


1.7 


-24 


.81 ± 


.18 


32. 


3 


0.33 


107 


3.8 ± 0.9 


x 


10 10 


00 


















NGC3923 


22.9 ± 


3.0 


-25 


.30 ± 


.28 


43. 


.8 


0.28 


47 


5.8 ± 1.6 


x 


10 10 


< 100 


L 


0. 


36 ± 


01 


47.9 ± 


0. 


5 


u 


NGC4125 


23.9 ± 


2.7 


-25 


.03 ± 


.25 


33. 





0.37 


80 


5.5 ± 1.4 


x 


10 10 


185 




0. 


.45 ± 0. 


01 


82.0 ± 


0. 
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a Distances and associated errors in Mpc, taken from Tonry ct al. (2001) and LEDA (Paturcl ct 
taken from the 2MASS extended source catalog dJarrett et al.j|2000l: absolute magnitude > 


al.ll997D, LEDA errors are assumed to be 15%. b Data 
(i-C-band), half-light radius in effective cllipticity 



and position angle (all J-band). c Absolute b lue luminosity in un its of Lb,q — 5.2 X 10 32 erg s ld Rotatio nal velocities in kms 1 evalu ated 
between 0.6 - 0.9 flj , t aken from Hyperled a (IPrugniel et al.lll998|) re fere nces. A: | Simien fc Prugniell U997al). B: IPrugniel fc Simienl dl994D . C : 

IBonfanti et al.l U995D , D:IBender et al.l Jl993) . E : ID Onofrio et al I U995D. F:IGraham et al.l 1119981). G: | Longo et ail J1994I) . H:ILonghetti et al.Nl998D. 
LICarollo fc Danziger] Ill994bn. J:IBettonil IU992T). K:IFisherl lU997fi . L: IPellegrini et al.l 1119971). M: | pavies fe Birkinshavd (I1988T) . N:| Simien fe Prugniell 
Il99sh . O: ISimien fe Prugniell lU997bT ). P: IBinnev et al.l dl99oK 7 O: ICarollo fc Danzigerl lll994ah . R: ICarollo et al.l fl993ll. S: "simien fe Prugniell 
(2002). ° Ellipticity an d position angle, eval uated between 0.6 — 0.9 Rj from published optical surface photometry. T: IGoudfrooii et all (I1994D ■ U: 

IBender et"afl U988I) . V: IPeletier et al.l U990D , W: IFalcdn-Barroso et~aLl ( 120061 ). 
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Fig. 2. — Adaptively binned Chandra X-ray gas surface brightness maps (right) and optical DSS ij-band images (left). The objects are 
ordered by 2MASS K band luminosity, starting with the most luminous galaxy on the top left of the page, and decreasing to the right and 
then to the next row. The physical scale of each image is 50 kpc X 50 kpc, except where indicated by an individually attached scale bar. 
The color range of the X-ray gas distribution is scaled logarithmically between 5 X 10 -11 and 3 X 10 — 7 photons sec - 'cm~ 2 arcsec - 2 . The 
x and y-axes are labelled according to right ascension and declination (2000), respectively. 




Fig. 2. — Continued 



The Chandra Gas Gallery of Normal Ellipticals 



11 




Fig. 2. — Continued 



12 



Diehl & Statler 




2 25:012:24:5B 



Fig. 3. — Adaptively binned Chandra X-ray gas surface bright- 
ness maps (right) and combined 2MASS J, H and _ft" s -band images 
(left). Figures [3^-ag showing all 33 galaxies with sufficient data to 
compute the X-ray ellipticity ex are available in high-resolution in 
the electronic edition of the Journal. The printed edition contains 
only a sample (NGC 4374). When available, 2MASS images were 
taken from the 2MASS Large Galaxy Atlas llJarrett et al.l 120031 . 
http : //irsa. ipac . caltech. edu/applications/2MASS/LGA/at las .html 
otherwise ima ges are taken f rom the 2MASS Extended Source 
Catalog (Jarrett ct al. 2000). The physical scale of each 
image is 6i?j across. The color range of the X-ray gas dis- 
tribution is scaled logarithmically between 5 X 10 — 1 and 



1 X 10 "photons sec cm 



and is identical to that 



of Figure [2] The x and y-axes are labelled according to right 
ascension and declination (2000), respectively. 
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Fig. 4. — Total X-ray gas luminosity in ergs as a function of 
absolute blue luminosity in units of blue solar luminosity (Lg q = 
5.2 X lO^ ergs -1 ) taken from LE DA. The dashed line shows the 
best fit by O'Sullivan et al. (2001), which they statistically correct 
for the expected contribution from unresolved point sources. Our 
data are consistent with their relation. 

photometry, and are not plotted in Figure removing 
the remaining two (NGC 4374 and NGC 5044) makes 
an insignificant difference. Excluding group member 
galaxies, or systems with significant optical isophotal 
twists or ellipticity gradients also has negligible effect. 

This lack of correlation is very surprising if hydrostatic 
equilibrium holds precisely. Roughly speaking, if one as- 
sumes an oblate logarithmic potential, the equipotential 
surfaces should be about one third as flattened as the un - 
der lying density distribution (|Binnev fc Trema inc 1987). 
Instead, we find that the X-ray isophotes are often much 
flatter than the stellar isophotes. We fully discuss the 
implications of this result, with detailed modeling, in § 
14.11 below. 

We notice a weak tendency for galaxies to have their 
X-ray gas major axis aligned with the stellar distribu- 
tion. Given the fact that X-ray isophotal twists often 



correlate with the orientation of asymmetries in the gas, 
it is not clear whether this alignment is a consequence of 
the underlying potential. An alternative is that it is a 
consequence of a misalignment between radio sources op- 
tical major axes (jPalimaka et al.lll979l ). which is present 
in our sample and which we fully discuss in Paper II. 

3.3. Rotational Support 

A significant fraction of the hot gas in elliptical 
galaxies is believed to be d ue to stellar mass-loss (e.g. 
iBrighenti fc Mathews) I2OOO0 . Thus, the gas should con- 
tain a significant amount of specific angular momen- 
tum, causing it to settle int o a flattened rotating disk 
(|Brighenti fc Mathews! Il997t ). This scenario predicts a 
relationship between the stellar rotational velocity and 
),the X-ray gas ellipticity, with faster rotating systems ex- 
hibiting flatter gas isophotes. Figure [7] shows the ob- 
served X-ray ellipticity between 0.6 and 0.9i?j as a func- 
tion of mean stellar rotational velocity in the same radial 
range. In agreement with previous ROSAT and Einstein 
obser vations of elliptical galaxies ijHanlan fc Bregmanl 
2000), we find no evidence for any correlation between 
these parameters. A Spearman rank analysis yields a 
probability of 54% for the null hypothesis of no corre- 
lation. While this is formally an upper bound on the 
probability due the presence of finite measurement errors 
and the presence of upper limits, it is obvious from the 
figure that the errors are not responsible for the absence 
of correlation. 

Such a lack of correlation could still be consistent with 
an alternative scenario, where the majority of the hot 
gas is not due to stellar mass loss, but instead is ac- 
quired externally through mergers, infalling gas clouds 
or through tidal stripping during close encounters. The 
hot gas would then slowly fl ow toward the center and set - 
tle into a cooling disk (e.g. Bri ghenti fc Mathews! 1 19971 ). 
Thus, if angular momentum is conserved during the in- 
flow, one should be able to observe a trend for the X -ray 
ellipticities to rise inward (Hanl an fc Bregmanll2000T ). as 
the angular momentum becomes more important. The 
detailed ellipticity profiles in Figure [5] do not support 
this prediction observationally, since ellipticities do not 
systematically increase inward. Moreover, significant ro- 
tational support would not account for the prevalence of 
strong isophotal twists. 

4. DISCUSSION 

4.1. Implications for Hydrostatic Equilibrium 

Within one stellar effective radius, the gravitational 
potentials of normal elliptical galaxies are dominated 
by the stellar component. Although some studies 
have found dark matter contributions in individual ob- 
jects as high as 40% within this radius (|Gerhard et al.l 
120011: iPadmanabhan et ail I2004D . in general the dark 
matter halo contribu tion should be secondary (e.g. 
iMamon fc Loka s 2005). Thus, if the assumption of hy- 
drostatic equilibrium were valid, one would expect the 
hot gas isophotes to reflect the shape of the stellar po- 
tential (|Buote fc Canizareill994h . Since projected isopo- 
tentials are rounder than the projected isopleths of the 
underlying density, we would also expect the gas emis- 
sion to be rounder than the starlight. However, Figure 
[S] shows that the gas is flatter than the starlight in one- 
third of the sample, and that gas and optical ellipticities 
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Fig. 5. — Radial isophotal ellipticity and position angle profiles for the subset of 36 elliptical galaxies containing sufficient signal. Large 
filled circles with error bars denote X-ray gas profiles; other symbols denote optical profiles as indicated in the figure. The 2MASS ellipticity 
and position angle arc marked with arrows at the right border. The vertical dashed lines mark the 2MASS J-band effective radius. 
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Fig. 5. — Continued. 
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Fig. 6. — X-ray gas ellipticity vs. optical cllipticity, both eval- 
uated between 0.6 — 0.9 Rj. There is no correlation, contrary to 
what what would be expected if the gas were in hydrostatic equi- 
librium in a stellar-mass-dominated potential. See also Figure [8] 
for theoretical predictions. 
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Fig. 7. — X-ray ellipticity as a function of rotational stellar 
velocity between 0.6 — 0.9 Rj. No correlation is present, contrary 
to what is expected if the hot gas would be due to stellar mass-loss 
and subsequently settle into a cooling disk. 

are completely uncorrelated. In fact, the Spearman rank- 
order correlation coefficient is slightly negative, with a 
value of —0.191. A two-sided confidence test shows that 
the difference from zero is not statistically significant. 10 
To determine whether the gas may still be hydro- 
static despite being morphologically uncorrelated with 
the starlight, we simulate hydrostatic systems in a va- 
riety of stellar-plus-dark-matter potentials. The stellar 
component is represented by a triaxial Hernquist model, 
whose mass density is given by 



p* = 



PH 



(4) 



(a/<z*)(l + a/a*) 3 ' 
where a is the semimajor axis of the ellipsoidal equiden- 

10 This statement also holds for extraction annuli ranging be- 
tween 0.5 and 1.2 as well as for the sample with the disturbed 
objects removed. 



sity surface, a* is the scale length, and pn is a normal- 
ization constant. The dark halo is taken to be either a 
triaxial NFW model, having density 



Pdm 



Pnfw 



(o/odm)(1 + a/auM) 2 



or a similar profile with a flat core: 



PDM 



Pnfw 



(1 + a/a D M) 



(5) 



(6) 



We parametrize the models by the short-to-long axis ra- 
tio c/a and triaxiality T = (a 2 — b 2 ) / (a 2 — c 2 ) of the stellar 
and dark matter components, the ratio of scale lengths 
a = (Zdm/a*) and the ratio M of dark to luminous mass 
inside the half-mass radius of the stellar component. 11 
The potenti als are computed u sing the formulae in the 
appendix of iFlores et al.l (|2005l ). 

The principal axes of the two components are assumed 
to be aligned and not intrinsically twisted. We do not 
consider misaligned or twisted models because there are 
no self-consistent equilibrium models for such configura- 
tions. Intrinsically misaligned or twisted galaxies, if they 
exist, would likely be in some normal mode of oscillation 
whose period is of order the dynamical time, and so any 
gas in these systems would not be hydrostatic. We also 
do not consider stellar distributions with radially vary- 
ing axis ratios. Even though this would tend to slightly 
decouple the local ellipticities of the starlight and the 
stellar potential, this is a moot point. Only two objects 
in Figure [5] show significant ellipticity gradients, and the 
results are unchanged if these are removed. 

The a priori axis ratio distribution of the stellar com- 
ponents is chosen to be approximately consistent with 
the observed ellipticity distribution of normal ellipticals. 
For populations of purely oblate or purely prolate sys- 
tems, we use rough analytic fits to the distributions ob- 
tained by I Vincent fe Rvdenl (|2005[ ) from the Sloa n Digi- 
tal Sky Survey Data Release 3. It is well known (|Rvdenl 
Il996t iTrcmbla v fe Merrittlll995T ) that purely axisymmet- 
ric populations cannot reproduce the observed ellipticity 
distribution without invoking negative numbers of galax- 
ies. We force our fits to be strictly positive, even though 
this means that the observed ellipticity distribution is 
not exactly recovered. For triaxial populations we use 
an analytic fit t o the "maximal ignorance" distribution of 
iBak fc StatTeil (|2000f ). This distribution is derived from 
simultaneous fits to surface photometry and stellar kine- 
matics for 13 galaxies, and is consistent with the overall 
ellipticity distribution. 

The combined potential is filled with isothermal gas, at 
a temperature that mimics the typical observed surface 
brightness profiles. We orient the system at random, and 
project the X-ray emissivity, which we take to be propor- 
tional to the square of the gas density. We calculate the 
stellar and gas isophotal ellipticities at 0.75 of the stel- 
lar effective radius, 12 to correspond with Figure [6l The 
results are quite insensitive to the assumed temperature; 

11 The half-mass radius of the Hernquist model is (l + 2 1 / 2 ) a*, 
the semimajor axis of the isodensity surface containing half the 
mass. 

12 The effective radius is the semimajor axis of the isophote 
enclosing half the projected light, which depends on the axis ratios 
and viewing geometry. 
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we are therefore confident that they will apply to mod- 
erately non-isothermal equations of state as well. 

We first compute models in which the shapes of the 
dark halos are correlated with the shapes of the stellar 
compone nts. The latter are ta ken to be triaxial, drawn 
from the iBak fe Statleri <|2000l ) distribution. For each 
object, the halo triaxiality T and axis ratio c/a are cho- 
sen randomly from a gaussian distribution centered on 
the stellar values, with a width a. We consider strongly 
(cr = 0.2) and weakly (a = 0.4) correlated halos. We 
compute an ensemble of several hundred models, oriented 
at random, to define the predicted distribution in the ex 
- e op t plane. Examples for strongly and weakly correlated 
halos are shown in Figure [5^ and [SJd. Not surprisingly, 
these models predict a significant correlation between ex 
and e op t. 

To assess quantitatively whether the models can be 
consistent with the observed distribution, we create sev- 
eral thousand simulated 24-object samples. For each sim- 
ulated sample, we randomly perturb the ex and e pt val- 
ues according to the actual measurement errors in the 
observed sample. We then compute the Spearman rank- 
order correlation coefficient and the unweighted mean 
ellipticity {ex), and determine the joint probability that 
the model will produce a Spearman coefficient as low as, 
and a mean ellipticity as high as, the observed values. In 
assessing the models we consider M values of 0.1, 0.3, 1, 
3, and 10, and a values of 1, 3, and 10. 

We find that all of the correlated halo models are ruled 
out at > 99.9% confidence. This is the case even if dark 
mass exceeds stellar mass by as much as a factor of 10 in- 
side the stellar half-mass radius. Clearly, if the isophotal 
ellipticity of hydrostatic gas is to be uncorrelated with 
that of the starlight, the shape of the halo must be un- 
correlated with that of the stars. 

Accordingly, we next consider a variety of uncorrelated 
halo models. These come in two types: (1) triaxial stellar 
components with either purely oblate or purely prolate 
halos having flattenings c/a of 0.5 and 0.2; and (2) ax- 
isymmetric (oblate or prolate) stellar components with 
random triaxial halos. Examples are shown in Figure 
[8};-f. Here we find that all of the following are excluded 
at > 99% confidence: (1) oblate halos; (2) random tri- 
axial halos with mass ratios M < 3; (3) prolate halos 
with c/a = 0.5_and M < 3; and (4) prolate halos with 
c/a = 0.2 and M < 0.3. 

This result does not rely on the presence of a large 
number of clearly disturbed objects in the sample. Six 
of the eight most asymmetric systems are not represented 
in Figure [6] because they do not have published optical 
surface photometry. Nonetheless, the objects with the 
largest ex values in Figure [5] are NGC 4374 and NGC 
5044. Of these, the former is clearly disturbed by its 
radio source, and the latter is lopsided at large radius, 
suggesting an interaction with the intergalactic medium. 
If we cull these two objects from the sample, the con- 
straints are weakened slightly but our qualitative result 
still holds. Correlated halo models are still ruled out 
at > 99.9% confidence, and random triaxial halos with 
mass ratios M < 3 are still excluded at 99% confidence. 
Oblate halos are excluded at 97% confidence, and pro- 
late halos with c/a = 0.5 and M < 3 are excluded at 
98% confidence. Prolate halos with c/a = 0.2 cannot be 
excluded for M > 0.1. 



Thus the only viable hydrostatic models are those in 
which the shape of the halo has nothing to do with the 
shape of the luminous galaxy; moreover, the dark halo 
must either dominate over luminous matter inside the 
stellar half-mass radius or be prominently sausage- or 
cigar-chaped. Both of these options conflict with well 
established results. There is widespread agreement, 
from stellar population modeling coupled with stellar 
dynamics or X-ray data, that normal ellipticals are 
not dar k-matter-dominated i n the i nner 1 effective 
radius (lLoewenstein fc White! Il999t iMamon fc Lokad 
l2005HLintott et al.ll2006HHumphrev fc Buotdl2006fl The 
mass fractions obtained in these studies correspond to 
0.2 < M < 1.0, which are ruled out by our simula- 
tions except for the case of fiat, prolate halos. But 
numerous s imulations of galaxy merg e rs and structure 
form ati on dBarnes fc Hernquistl 11996 : ISugerman et al.l 
20001 iMeza et all 120031: iBailin fc Steinmetzl 120051: 



Allgood et al.l 120061 : ICox et alJ l2006h demonstrate that 



halos should be neither prolate nor highly flattened. 
Moreover, there is no compelling reason to expect 
that dark and stellar matter distributions should have 
very different shapes; both obey the same collisionless 
dynamics, and should end up similarly distributed 
following a major event such as a merger. We are 
therefore forced to the conclusion that the assumption 
of hydrostatic equilibrium must be incorrect. 

To be more precise, even though the gas in normal el- 
lipticals may be, on average, close enough to hydrostatic 
that reasonable radial mass profiles can be derived, hy- 
drostatic equilibrium does not hold tightly enough that 
X-ray morphology can be used to trace the shape of 
the total mass di s tribution. Buote and collabo rators 
(|Buote et al.ll2002l : iBuote fc Canizareslll996l . fl99l have 
used this technique to argue for a flattened dark matter 
halo in NGC 720. While it is possible that hydrostatic 
equilibrium could hold in some individual systems, our 
results show that this would be the exception rather than 
the rule, and there is no independent test to reveal where 
hydrostatic equilibrium does hold. Also, in the case of 
NGC 720, we find that properly removing the resolved 
and unresolved LMXB emission reduces the gas ellip- 
ticity, and, along with it, the need for a flattened dark 
halo. The ex profile we show in Figure [5] appears to 
be entire ly consistent with the oblate mass-follows-light 
model of IBuote etakl (|2002l Fig. 3) at radii between 10" 
and 80". Ho wever, a direc t comp arison is complicated by 
the fact that IBuote et~a l. (2002) do not measure isopho- 
tal ellipticities and position angles, but rather their mean 
values inside elliptical apertures. 13 

To preserve hydrostatic equilibrium in the general 
population of ellipticals, our data would require ex- 
tremely flat, complex dark matter halos with signifi- 
cant substructure to explain the large ellipticities and 
isophotal twists. This is highly unlikely and inconsis- 
tent with stellar kinematics (e.g. ICappellari et al.ll2006f ) 
or simulations of l arge-scale structure formation (e.g. 
iSpringel et alJ 2005). Even though triaxial dark matter 
halos are able to produce simple t wists in X-ray isophotcs 
(|Romanowskv fc Kochaneklll998f ). these effects would be 

13 The mass-follows-light model does not reproduce the offset 
between optical and X-ray position angles, a failing that does not 
have a unique interpretation. 
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insufficient to explain the complexity and magnitude of 
the observed features. 

We see no delineation between disturbed systems that 
are not hydrostatic and undisturbed systems that are. 
While we cannot prove that the round, symmetric sys- 
tems are not hydrostatic, we emphasize that the sample, 
as a whole, has failed the test. If one chooses to invoke 
hydrostatic equilibrium for a system that appears "re- 
laxed," it should be clearly understood that this reflects 
a hope rather than a verifiable assumption. 

We take the view that the isophotal ellipticity of the 
X-ray gas does not reflect the underlying potential. El- 
lipticity is merely a measure of the m — 2 Fourier am- 
plitude of the deviations from circular symmetry. In the 
optical, the other Fourier harmonics are small by com- 
parison; in the X-ray, they are not. In Paper II of this 
series we will argue that the gas ellipticity is a measure 
of morphological disturbance, quantify that disturbance, 
and demonstrate that it is likely produced by AGN. 

The assumption of hydrostatic equilibr ium is widely 
used t o derive radial ma s s profi le s (e.g. iForman et al.l 
1981 iKilieen fc Bicknelj [1981 iPaolillo et all 1200? 
Humphrey et alJl200d iFukazawa et al.ll2006l ). If the gas 
is not hydrostatic, what are the consequences? Imagine 
a region that was locally overpressured by a factor q. As- 
suming adiabatic expansion this region would expand by 
a factor q 1 ^ in linear size over a sound crossing time scale 
of ~ 10 8 years, to regain pressure equilibrium. To redis- 
tribute gas a significant distance along an equipotential 
would thus imply an overpressure q < 10. If the galaxy 
were globally overpressured by this factor, one would in- 



fer a mass a factor of q too large. However, normal el- 
lipticals are probably only locally overpressured, and so 
averaging azimuthally over N over- and underpressured 
regions in a radial analysis would tend to reduce this 
effect by a factor of order N 1 / 2 , resulting in a mass es- 
timate in error by at most a factor of a few, either high 
or low. Subsonic bulk motions generated by such tem- 
porary overpressures could also enhance or reduce radial 
pressure support, de pending on their radial gradient. 

An early study bv lSchindlerl ([1996) demonstrated that 
mass profiles of simulated clusters can be recovered suc- 
cessfully with the hydrostatic equilibrium method, lim- 
ited primarily by the presence of significant substruc- 
ture. More recent studies of hydrodynamic cluster for- 
mation simulations (jRasia et al.112004 : iNagai et al.ll2007T ) 
show that neglecting bulk motion in a hydrostatic anal- 
ysis t ends to underestim ate the gravitating mass by 5 to 
20%. iRasia et"aTI (2006) find a further systematic effect 
of similar magnitude that they attribute to the complex 
temperature structure of the cluster; however this ap- 
pears only when they reduce the backgrounds in their 
simulated observations far below typical Chandra values. 
Unfortunately, no such study is available for the mass 
regime of normal elliptical galaxies and none yet takes 
into account disturbances due to AGN activity, which our 
simple arguments above suggest may be larger in magni - 
tude. On the observational side, iFukazawa et al.1 (2006) 
derive total mass profiles from Chandra X-ray data of 53 
normal elliptical galaxies assuming hydrostatic equilib- 
rium and compare 7 of these to mass profiles based on 
stellar kinematics. They find that for 2 out of 7 galaxies, 
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total enclosed masses (inside observationally accessible 
maximum radii) differ by a factor of ~ 2 and for one 
object, NGC 3379, the diff e rence even grows to a factor 
of - 7. iPellegrini fc Croft] (|200ll reconcile this discrep- 
ancy for NGC 3379 by fitting its X-ray properties with a 
wind model, staying consistent with the optically derived 
mass profile. They conclude that the gas is not in hy- 
drostatic equilibrium, but in a general state of outflow, 
calling the X-ray mass measurement into question. 

4.2. Implications for Rotational Support 

The hot gas in elliptical galaxies is believed to come 
from some combination of stellar mass loss, infall or 
mergers. In each case, the gas should carry sig- 
nificant amounts of angular momentum. In a stan- 
dard cooling flow model that conserves angular mo- 
mentum, the gas will settle into rotationally sup- 
ported cooling disks (iBrighenti fc Mathews! 119971 119961 : 
Kiev fc Mathewsl 119951) . These rotating cooling flow 
models predict cllipticity profiles that rise inward, with 
minimal isophotal twist. 

Figures [5] and [7] show that these predictions are incon- 
sistent with observations. Although we do find rather 
large X-ray ellipticities, they do not systematically in- 
crease toward smaller radii. Our profiles also reveal dra- 
matic isophotal twists and asymmetries that would not 
arise naturally in a scenario where the ellipticities are 
caused by rotation alone. In addition, we find no rela- 
tion between stellar rotational velocity and the gas flat- 
tening, as one would expect if the gas is mainly due to 
stellar mass loss. Thus, we conclude that the X-ray gas 
morphology is not dictated by rot ation, in agreement 
with previous ROSAT observations (jHanlan fc Bregmanl 
120001) . 

The question now is whether the failure to observe ro- 
tationally flattened X-ray disks is a serious blow to cool- 
ing flow models. One way to save the cooling flow pic- 
ture would be to efficiently tran sfer angular momentum 
t hrou gh turbulence in the gas (IShadmehri fc Ghanbarl 
I2002t ferighenti fc Mathewsl |2000D . but it is unclear how 
effective this process is for elliptical galaxies. An- 
other alternative is that spatially distributed multi- 
temperatu re mass dropout can cir cularize the X-ray 
isophotes (Ha nlan fc Bregm an 2000). However, recent 
XM M-Newton spectroscopy of norma l elliptical galaxies 
(e.g lXu et al.ll200"2l : iBuote et al.i r2003) rules out the exis- 
tence o f sufficient intermedia te temperature gas. In ad- 
dition, iBregman et alJ (|2006l ) are able to directly trace 
the amount of this intermediate gas by measuring its 
O VI emission (~ 10 5 5 K) w jth the F ar Ultraviolet Spec- 
tral Explorer (FUSE). lBTegman et all (|2005l ) examine the 
spatial distribution of this warm gas for NGC 4636 and 
NGC 5846 and constrain the size of the cooling region to 
be smaller than 0.8 kpc and 0.5 kpc, respectively. Their 
measurements are consistent with moderate cooling flows 
with centrally concentrated mass dropout and rule out 
spatially distributed dropout for these galaxies. 

Another explanation for the lack of disk signatures is 
that the hot gas is actually flowing outward instead of 
inward, eliminating the ne ed for the gas to settle into a 
cooling disk. An idea by iBrighenti fc Mathewsl (|2006l ) 
seems to be successful in stopping inward cooling flows 
by AGN-induced massive jet outflows, which are stable 
for several gigayears. 



Alternatively, morphological asymmetries, which are 
clearly significant, may be masking possible disk signa- 
tures. Whatever the cause of these asymmetries (see Pa- 
per II), this effect could be effective in disrupting the 
flattest X-ray isophotes at small radii close to the center. 
Detailed simulations are needed to confirm this possi- 
bility. Thus, our results do not imply that rotational 
support of the hot gas in elliptical galaxies has no im- 
portance at all. We can merely state that it is not the 
dominant factor that causes the large observed elliptici- 
ties. 

There is other evid ence that gas rotatio n migh t be 
present. For example, IStatler fc McNamaral (|2002f ) find 
extremely flat X-ray emission in NGC 1700, which they 
interpret as a large-scale cooling disk. Their model yields 
a specific angular momentum and cooling time for the 
hot gas that is consistent with the gas having been ac- 
quired during the last major merger. We find that re- 
moval of unresolved point sources somewhat reduces the 
X-ray ellipticity of this object; but the data are still con- 
sistent with their claimed 15 kpc rotating disk. Addi- 
tional support for rot ation has recently been reported by 
iBre gman et all (|2005[ ). who find that the O VI line struc- 
ture of NGC 4636 exhibits signs consistent with rotation. 

5. CONCLUSIONS 

We have analyzed the X-ray emission of 54 normal el- 
liptical galaxies in the Chandra archive and isolated their 
diffuse hot interstellar gas emission from the emission of 
discrete stellar point sources for the first time. We qual- 
itatively and quantitatively compare the morphology of 
the hot gas to the shape of the stellar distribution and 
find that they have very little in common, despite the 
known Lx~Lb relation. We compute ellipticity and po- 
sition angle profiles for the X-ray gas and compare them 
to published optical profiles. In particular, we do not find 
a correlation between optical and X-ray ellipticities mea- 
sured in the regions where stellar mass dominates the po- 
tential, suggesting that the gas is at least far enough out 
of hydrostatic equilibrium that the information about the 
shape of the underlying potential is lost. We also argue 
that X-ray derived radial mass profiles may be in error 
by factors of as much as a few, without the necessity for 
the galaxy to drive a global wind. 

Although we find large X-ray gas ellipticities to be 
common, the gas morphology is generally inconsistent 
with rotationally flattened disks, and a comparison with 
stellar rotational velocities yields no evidence for signifi- 
cant rotational support. 

The fact that neither the shape of the underlying po- 
tential nor rotational support determine the overall dis- 
tribution of the X-ray emitting gas, combined with its 
general disturbed appearance, suggests the involvement 
of another major component: the central AGN. We as- 
sess the importance of the AGN in Paper II, where we 
draw a connection between gas morphology and AGN 
luminosity. These new findings arc consistent with the 
AGN constantly stirring up the interstellar medium by 
inflating buoyant bubbles, which may also play a role 
in redistributing the angular momentum of the hot gas 
through entrainment. These intermittent AGN outbursts 
could also be responsible for disrupting or masking the 
signatures of cooling disks in the central regions of ellip- 
tical galaxies. 
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